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Abstract 

We have calculated spectral functions associated with hadronic current correlation functions for 
vector currents at finite temperature. We made use of a model with chiral symmetry, temperature- 
dependent coupling constants and temperature-dependent momentum cutoff parameters. Our 
model has two parameters which are used to fix the magnitude and position of the large peak seen 
in the spectral functions. In our earlier work, good fits were obtained for the spectral functions 
that were extracted from lattice data by means of the maximum entropy method (MEM). In the 
present work we extend our calculations and provide values for the three-momentum dependence 
of the vector correlation function at T = 1.5 T c . These results are used to obtain the correlation 
function in coordinate space, which is usually parametrized in terms of a screening mass. Our 
results for the three-momentum dependence of the spectral functions are similar to those found 
in a recent lattice QCD calculation for charmonium [S. Datta, F. Karsch, P. Petreczky and I. 
Wetzorke, hep-lat/0312037 . For a limited range we find the exponential behavior in coordinate 
space that is usually obtained for the spectral function for T > T c and which allows for the definition 
of a screening mass. 
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I. INTRODUCTION 



In a number of recent works [1- 3] we have calculated various hadronic correlation func- 
tions and compared our results to results obtained in lattice simulations of QCD [4-6]. The 
lattice results for the correlators, G(t, T), may be used to obtain the corresponding spectral 
functions, a(u,T), by making use of the relation 

poo 

G(t,T)= dua( y u,T)K( y r,u,T) , (1.1) 
Jo 

where 

cosier -1/2T)] 
K ( T > U > T > = sinh(c/2T) • (L2) 

The procedure to obtain a(ou, T) from the knowledge of G(t, T) makes use of the maximum 
entropy method (MEM) [7-9], since G(r,T) is only known at a limited number of points. 

In our studies of meson spectra at T = and at T < T c we have made use of the Nambu- 
Jona-Lasinio (NJL) model. The Lagrangian of the generalized NJL model we have used in 
our studies is 



C = q(ij - m°)q + ^ ]T[(gA*g) 2 + (qt^q) 2 ] (1.3) 

i=0 

~ ^E[(^7^) 2 + (^ 7 57^) 2 ] 

i=0 

G 

+ _^{det[g(l + X 5 )q] + det[g(l - X 5 )q}} + C conf . 

Here, m° is a current quark mass matrix, m°=diag(m°, m°, m°). The \ are the Gell- 
Mann (flavor) matrices and A = a/2/31, with 1 being the unit matrix. The fourth term is 
the 't Hooft interaction and C con f represents the model of confinement used in our studies 
of meson properties. 

In the study of hadronic current correlators it is important to use a model which respects 
chiral symmetry, when m° = 0. Therefore, we make use of the Lagrangian of Eq. (1.3), while 
neglecting the 't Hooft interaction and C con f. Thus, there are essentially three parameters 
to consider, Gs, Gy and a cutoff parameter k max , which restricts the momentum integrals 
so that k < k max . When we use the NJL model to study matter at finite temperature, we 
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introduce the temperature-dependent parameters Gs(T) and GV(T). (We have also used a 
Gaussian cutoff for the momentum integrals in our earlier work.) These parameters have 
been adjusted to obtain fits to the spectral functions a(uj, T) for T/T c = 1.5 and T/T c = 3.0, 
which are the values of T/T c studied in the lattice simulations of QCD [10]. 

The temperature-dependent coupling constants and cutoff parameters of our work are 
analogous to the corresponding density-dependent parameters introduced in Ref. [11] and 
[12]. Further study of models with temperature-dependent and density-dependent parame- 
ters are of interest and a general theoretical formalism for the introduction of such depen- 
dencies should be considered. 

In this work we limit our study to the data for the vector correlator at T = 1.5 T c [10, 
13] and therefore only need to specify Gy(T) and the momentum cutoff. (We remark that 
the results for the scalar, vector, pseudoscalar and axialvector correlators are quite similar 
in the deconfined region.) In Figs. 1 and 2 we show the data obtained by the MEM method 
at T/T c = 1.5 and T/T c = 3.0 for both pseudoscalar and vector correlators [10, 13]. The 
second peaks in these correlators are known to be a lattice artifacts [13]. 

The organization of our work is as follows. In Appendix A we present the formalism for 
calculation of pseudoscalar and vector correlation functions. In Appendix A we discuss the 
calculation of the correlator in the case the quark and antiquark carry zero total momentum. 
In Appendix B we show how the formalism is modified for the correlator calculated at finite 
momentum P . In Section II we present the results of our calculation of the imaginary part 
of the correlator a(cu, P). (Since we place P along the z-axis this quantity may be written 
as a(uj,0,0, P z ) in accord with the notation of Ref. [10].) Our results for a(u, P) will be 
presented for a series of values of | P | in Section II. In Section III we present our result for 
the coordinate-dependent correlator C(z) which is proportional to the correlator defined in 



Eq. (1) of Ref. [10], 




(1.4) 



We may also use the form 



C(z) 



4 



1 



— oo 



r 

dP z e lP * z 



Jo 



oo 



dp 



, 2 a(P 2 ,0,0,P z ) 

p2 



(1.5) 



Finally, in Section IV we present our conclusion and further discussion. 
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FIG. 1: The spectral functions a/uj 2 for pseudoscalar states obtained by MEM are shown [10, 13]. 
The solid line is for T/T c = 1.5 and the dashed line is for T/T c = 3.0. The second peak is a lattice 
artifact [13]. 

II. TEMPERATURE DEPENDENT HADRONIC CURRENT CORRELATORS 
AT FINITE MOMENTUM 

We make use of the formalism presented in Appendices A and B to obtain values of the 
vector correlator at T = 1.5 T c . Here we take T c = 270 MeV, since we have usually made 
comparison to lattice calculations made in the quenched approximation. 

In Fig. 3 we present <j(uj)/uo 2 , for various values of |P|, as function of uo 2 . Comparison 
may be made to the lattice data shown in Fig. 2 [10, 13]. (We again note that our calculation 
does not reproduce the second peak in the lattice data which is known to be a lattice artifact 
[13].) The curves shown in Fig. 3 are given for values of | P | ranging from 0.10 GeV to 2.10 
GeV in steps of 0.20 GeV. For these calculations we have used k max = 1.22 GeV. Our results 
for the various values of | P \ given in Fig. 3 may be compared to Fig. 20 of Ref . [14] . We see 
that the results calculated by completely different methods are similar. 
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FIG. 2: The spectral functions a/uj 2 for vector states obtained by MEM are shown [10, 13]. See 
the caption of Fig. 1. The second peak is a lattice artifact [13]. 



0.20-, 
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FIG. 3: The imaginary part of the correlator a(uo)/uj 2 is shown for various values of |P| as a 
function of lo 2 . Starting with the topmost curve the values of \P\ in GeV units are 0.10, 0.30, 
0.50, 0.70, 0.90, 1.10, 1.30, 1.50, 1.70, 1.90 and 2.10. Here we have used G s = 1.2 GeV~ 2 and 
k max = 1.22 GeV. 
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FIG. 4: The correlation function C(z) denned in Eq. (1.5) is shown. The dotted line represents a 
fit using an exponential function. 

III. THE COORDINATE-SPACE CORRELATION FUNCTION 

We have used the results of our calculations which were presented in Fig. 3 to calculate 
C(z) of Eq. (1.5). The result of that calculation is shown in Fig. 4. We note that the simple 
assumption for the behavior of this correlator that is usually made, C(z) ~ exp[— m sc z], is 
born out in our calculation for 1 GeV -1 < z < 3 GeV -1 . For the study of charmonium on 
the lattice, values found for the screening mass are given in Ref. [14]. 

For the result shown in Fig. 4 we obtain a screening mass of 1.02 GeV which may be 
compared to the value of 7rT=1.27 GeV. The dotted line represents an exponential fit to our 
result. 

IV. DISCUSSION 

Recent theoretical work concerning the quark-gluon plasma has been discussed by 
Shuryak Q|. He notes that the physics of excited matter produced in heavy- ion colli- 
sions in the region T c < T < 4T C is different from that of a weakly coupled quark-gluon 
plasma. That is due to the strong coupling generated by the bound states of the quasipar- 
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tides. These bound states appear as resonant structures when the imaginary parts of the 
hadronic correlators are extracted from lattice data using the maximum entropy method 
(MEM) [4-9]. These resonances lead to very strong interactions between the quasiparticles 
and are, in part, responsible for very small mean free paths and collective flow in heavy- 
ion collisions. That flow may be described by hydrodynamics. Indeed, Shuryak suggests 
that "... if the system is macroscopically large, then its description via thermodynamics of 
its bulk properties (like matter composition) and hydrodynamics for space-time evolution 
should work" 

In our present work we have presented a simple chiral model that is able to reproduce the 
resonances extracted from the lattice data via the MEM procedure. We have calculated the 
imaginary parts of the vector correlator for various values of the total external momentum 
P. We have also calculated the correlation function, C(z), and find the simple exponential 
behavior C(z) ~ exp[— m sc z] for a limited range of z. We have shown in an earlier work that 
exponential behavior with the appropriate screening mass may be obtained for the full range 



of z values if quite small values of k max (of the order of 0.4 GeV) are used [2Q]. Finally, we 
note that an extensive discussion of screening masses appears in Ref. 



APPENDIX A 

For ease of reference, we present a discussion of our calculation of hadronic current 
correlators taken from Ref. [3]. The procedure we adopt is based upon the real-time finite- 
temperature formalism, in which the imaginary part of the polarization function may be 
calculated. Then, the real part of the function is obtained using a dispersion relation. The 
result we need for this work has been already given in the work of Kobes and Semenoff 



17j . (In Ref. the quark momentum is k and the antiquark momentum is k — P. We 



wil 



17 



adopt that notation in this section for ease of reference to the results presented in Ref. 



.) With reference to Eq. (5.4) of Ref. 17|. we write the imaginary part of the scalar 
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polarization function as 

Im J S {P\ T) = -N C (3 S e(P°) / „t ttt (Al) 

M ' ; 2 MA v V (2tt) 3 \2E x {k)2E 2 {k) J v ; 

x{[l - m(A;) - n 2 (A;)]5(P - ^(fc) - £ 2 (fc)) 
-M*) - n 2 (k)]5(P° + Erik) - E 2 {k)) 
~[n 2 (k) - ni (k)}5(P° - Ex(k) + E 2 (k)) 
-[1 - m(fc) - n 2 (k)}5(P° + E x {k) + • 

Here, E x {k) = [k 2 + m\(T)} 1 / 2 . Relative to Eq. (5.4) of Ref. jl7|. we have changed the sign, 
removed a factor of g 2 and have included a statistical factor of N c . In addition, we have 
used a sharp cutoff, k max , for the momentum integral. We also note that 



= eP E l(k) + V ( A2 ) 



and 



H2 ^ = e PE 2 (k) + x • ( A3 ) 

For the calculation of the imaginary part of the polarization function, we may put k 2 = 
m\(T) and {k — P) 2 = m 2 (T), since in that calculation the quark and antiquark are on- 
mass-shell. In Eq. (Al) the factor (3g arises from a trace involving Dirac matrices, such 
that 

(3 S = -Tt[($ + mi )(| - f + m 2 )\ (A4) 
= 2P 2 - 2(mi + m 2 f , (A5) 

where m\ and m 2 depend upon temperature. In the frame where P = 0, and in the case 
mi = Tn 2 , we have (3s = 2P 2 (1 — 4m 2 /Pq). For the scalar case, with mi = m 2 , we find 

Im J S {P\T) = ^ (l - ^i) 3/2 [1 - 2m(A0] , (A6) 



where 



4 

with fc < 



fc 2 = ^l- m 2 (T), (A7) 
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For pseudoscalar mesons, we replace fts by 

P P = -Tr[i l5 (^ + mi )i l5 (lt-f + m 2 )} (A8) 
= 2P 2 - 2(mi - m 2 f , (A9) 

which for m x = m 2 is ftp = 2P 2 in the frame where P = 0. We find, for the ir mesons, 

Im J P (P 2 ,T) = ^ (l " A ~^) 112 [1 " ^{k)\ , (MO) 

where k 2 = Pq/4: — m 2 (T), as above, with k < k max . Thus, we see that the phase space 
factor has an exponent of 1/2 corresponding to a s-wave amplitude. For the scalars, the 
exponent of the phase-space factor is 3/2, as seen in Eq. (A6). 
For a study of vector mesons we consider 

= Tr[ 7/ # + -f + m 2 )\ , (All) 

and calculate 

gri%, = 4[P 2 - ml - mj + 4 mi m 2 ] , (A12) 

which, in the equal- mass case, is equal to 4P 2 + 8m 2 (T), when mi = m 2 and P — 0. This 
result will be needed when we calculate the correlator of vector currents. Note that, for 
the elevated temperatures considered in this work, m u (T) = m d (T) is quite small, so that 
4P 2 + 8m 2 (T) can be approximated by 4P 2 , when we consider the vector current correlation 
functions. In that case, we have 

Im J v (P 2 , T) ~ hm J P (P 2 , T) . (A13) 

At this point it is useful to define functions that are not for k > k max : 

Im J P (P 2 , T) = ^ ^1 - ) [1 - 2m(fc)] , (A14) 



N C P 2 ( 4m 2 (T)\ 1/2 

uiLjpyr- ,i) = 

and 

Im J V (P 2 , T) = 2 -^(l- V2 [1 - 2n 1 (Jb)] , (A15) 

For the functions defined in Eq. (A14) and (A15) we need to use a twice-subtracted dispersion 
relation to obtain Re J P (P 2 ,T), or Re J V (P 2 ,T). For example, 

P 2 

Re Jp(P 2 , T) = Re J P (0, T) + [Re J P (P 2 , T) - Re J P (0, T)] (A16) 
P 2 (P 2 -P 2 ) / A2 ^ Im Jp(s,T) 



K Jim^T) S(P 2 - S)(P 2 - S) ' 
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where A 2 can be quite large, since the integral over the imaginary part of the polarization 
function is now convergent. We may introduce Jp(P 2 ,T) and Jy(P 2 ,T) as complex func- 
tions, since we now have both the real and imaginary parts of these functions. We note that 
the construction of either Re Jp(P 2 ,T), or Re Jy(P 2 ,T), by means of a dispersion relation 
does not require a subtraction. We use these functions to define the complex functions 
J P (P 2 ,T) and J V (P 2 ,T). 

In order to make use of Eq. (A16), we need to specify Jp(0) and J P (P 2 ). We found it 
useful to take P 2 = -1.0 GeV 2 and to put J P (0) = J P (0) and J P (P 2 ) = Jp(P 2 ). The 
quantities Jy(0) and Jy(P 2 ) are determined in an analogous function. This procedure in 
which we fix the behavior of a function such as ReJy(P 2 ) or ReJy(P 2 ) is quite analogous 
to the procedure used in Ref. Q|. In that work we made use of dispersion relations to 
construct a continuous vector-isovector current correlation function which had the correct 
perturbative behavior for large P 2 — > — oo and also described the low-energy resonance 
present in the correlator due to the excitation of the p meson. In Ref. 18] the NJL model 
was shown to provide a quite satisfactory description of the low-energy resonant behavior 
of the vector-isovector correlation function. 

We now consider the calculation of temperature-dependent hadronic current correlation 
functions. The general form of the correlator is a transform of a time-ordered product of 
currents, 

iC(P 2 ,T) = J d 4 xe lP - x « T(j{x)j{0)) » , (A17) 

where the double bracket is a reminder that we are considering the finite temperature case. 

For the study of pseudoscalar states, we may consider currents of the form jp t i(x) = 
g(x)i75A*g(x), where, in the case of the n mesons, i = 1, 2 and 3. For the study of scalar- 
isoscalar mesons, we introduce js,i{x) = d(x)X l a(x), where i = for the flavor-singlet current 
and i = 8 for the flavor-octet current 1!) . 

In the case of the pseudoscalar-isovector mesons, the correlator may be expressed in terms 
of the basic vacuum polarization function of the NJL model, J P (P 2 ,T). Thus, 

Cr(P 2 .T) = MP\T) 1 _ Gp{T 1 )M ^ Ty (MS) 

where Gp(T) is the coupling constant appropriate for our study of 7r mesons. We have 
found Gp{T) = 13.49 GeV~ 2 by fitting the pion mass in a calculation made at T = 0, with 
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m u = rrid = 0.364 GeV. The result given in Eq. (A18) is only expected to be useful for small 
P 2 , since the Gaussian regulator strongly modifies the large P 2 behavior. Therefore, we 
suggest that the following form is useful, if we are to consider the larger values of P 2 . 



C P (P 2 ,T) 
P 2 



Jp{P\T) 
P 2 



l-G P (T)Jp(P 2 ,T)- 



(A19) 



(As usual, we put P = 0.) This form has two important features. At large Pq, 
Im Cp(P , T)/P 2 is a constant, since Im J p (Pq,T) is proportional to P 2 . Further, the 
denominator of Eq. (A19) goes to 1 for large Pq. On the other hand, at small P 2 , the 
denominator is capable of describing resonant enhancement of the correlation function. As 
we have seen, the results obtained when Eq. (A19) is used appear quite satisfactory. ( We 
may again refer to Ref. [IS], in which a similar approximation is described.) 

For a study of the vector-isovector correlators, we introduce conserved vector currents 
3n,i( x ) — <l( x )l^i<l(x) with i=l, 2 and 3. In this case we define 



JP(P 2 ,T)=[ 9 ^- 



P 2 



and 



CP(P 2 ,T)=lg 



piipu 

—p-2— 



Jv{P\T) 



C V (P\T) 



(A20) 



(A21) 



taking into account the fact that the current j^i(x) is conserved. We may then use the fact 
that 



J V (P 2 ,T) = ^ W J^(P 2 ,T) 



and 



lmJ v (P 2 ,T) 



P 2 + 2m 2 t (T) 
8vr 



-ImJ P (P 2 ,T). 



1 - 



4m 2 (T) 



p2 




1/2 



[i-2m(A0] 



(See Eq. (A7) for the specification of k — \k\.) We then have 

C V (P 2 ,T) = J V (P 2 ,T) 1 
where we have introduced 
ImJ y (P 2 ,T) 



P 2 + 2m 2 (T) 



l-G v (T)J v (P 2 ,T) 
4m 2 (T)\ 1/2 



7T 
2 



p2 
Ml 



7 



[i-2m{k)} 



lmJ P (P 2 ,T). 



(A22) 

(A23) 
(A24) 

(A25) 

(A26) 
(A27) 
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In the literature, ui is used instead of P Q [4-6] . We may define the spectral functions 

a v (u, T) = — Im C v (u, T) , (A28) 

7T 

and 

a P (uo, T) = - Im C P (u, T) , (A29) 

7T 

Since different conventions are used in the literature [4-6], we may use the notation 
ap(u, T) and av(u;, T) for the spectral functions given there. We have the following relations: 

a P {u,T) = a P (u,T), (A30) 

and 

a v (u,T) 3 

2 = 4^v(w, , (A31) 

where the factor 3/4 arises because, in Refs. [4-6], there is a division by 4, while we have 
divided by 3, as in Eq. (A22). 



APPENDIX B 

— * 

Here we extend the work of Appendix A to consider case of finite three-momentum, P. 
We consider the calculation of lmJ P (P°, P, T). The momenta P° and P are the values 
external to the loop diagram. Internal to the diagram, we have a quark of momentum 
k + P/2 leaving the left-hand vertex and an antiquark of momentum k — P/2 entering the 
left-hand vertex. It is useful to define 



Ei(k) = k + P/2 



p2 X 1/2 

= ( k 2 + — + kPcosOJ 



(Bl) 
(B2) 



and 



E 2 (k) = k - P/2 



= [k 2 + — - kPcosO 



1/2 



(B3) 
(B4) 



Here k 



and P = \P\ 
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We have 



lmJ v (P°,P,T) = -N c (3 v e(P° 



d 3 k 



2tt 



(2tt) 3 \2E 1 (k)2E 2 (k) 
x{[l - m(A;) - n 2 (k)]5(P° - E 1 (k) - E 2 (k)) 

-M*0 - n 2 (A;)]5(P° + Ei(k) - £ 2 (£;)) 

-M*0 - n^kMP - E ± (k) + £ 2 (£;)) 

-[1 - ni (k) - n 2 (k)]5(P° + E 1 (k) + £ 2 (*0)} . 



(B5) 



Here, 



and 



ni(k) = 



n 2 {k) 



(B6) 



(B7) 



e /3£ 2 (fc) _|_ 1 

In Eq. (B5), the second and third terms cancel and the fourth term does not contribute. It 
is useful to rewrite S(P° — E 1 (k) — E 2 {k)) using 

2 



where 



<S[/(cos0)] 



df 



-5(cos9 — x) , 



d cos 9 I x 



(B8) 



We find 



and obtain 



x 2 = 


-- cos 2 9 






AP 2 (k 2 + P 2 /4)-P 4 






4k 2 P 2 


df 


2 


E 1 (k)-E 2 (k) 


dcos9 


E\(k)E 2 (k) 



1 fkmax h-lAh. 

Im J P (P°,P,T) = l -N c p P e{P«)^) 2 J ^ 



I 



2E 1 (k)E 2 (k) 



[1 - ni(fc) - n 2 {k)\ 



df(cos9) 



dcos 9 
x5(cos9 — x)d(cos9) 



(B9) 



(BIO) 



(Bll) 
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We note there is a singularity when E\(k) = E2(k). That occurs when cos# = or 9 = tt/2. 
For our calculations we eliminate the point with 9 = 7r/2 when evaluating the angular 
integral over <i(cos6')5(cos6' — x) in the last expression. We obtain 



Im Jp(P°, P, T) = N c f3 P e(P ! 



Air 2 



(2tt) 3 
1 - rii(k) - n 2 (k) 



k 2 dk 



kP\E x {k) - E 2 (k)\ 



where x is obtained from Eq. (B9), 



x 



pO 

kF 



p2 p2 



1/2 



(B12) 



(B13) 
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